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ABSTRACT 

Many models of early structure formation predict a period of heating immediately preceding 
reionization, when X-rays raise the gas temperature above that of the cosmic microwave back- 
ground. These X-rays are often assumed to heat the intergalactic medium (IGM) uniformly, 
but in reality will heat the gas more strongly closer to the sources. We develop a framework 
for calculating fluctuations in the 21 cm brightness temperature that originate from this spatial 
variation in the heating rate. High-redshift sources are highly clustered, leading to significant 
gas temperature fluctuations (with fractional variations ^ 40%, peaking onk ^ O.lMpc"^ 
scales). This induces a distinctive peak-trough structure in the angle-averaged 21 cm power 
spectrum, which may be accessible to the proposed Square Kilometre Array. This signal 
reaches the ~ 10 mK level, and is stronger than that induced by Lya flux fluctuations. As well 
as probing the thermal evolution of the IGM before reionization, this 21 cm signal contains 
information about the spectra of the first X-ray sources. Finally, we consider disentangling 
temperature, density, and Lya flux fluctuations as functions of redshift. 
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1 INTRODUCTION 

The formation of the first luminous objects ends the cosmic 
"dark ages" and begins a period of heating and ionization of 
the intergalactic medium (IGM). The global thermodynamic 
history of this epoch, which culminates in reionization, de- 
pends upon many poorly constrained processes such as star 
formation, radiat i ve fe edback, and the growth of HII regions 
tearkana & Loebl 1200 ih . Currently, the best constraints on the 
ionization history come from obser vations of the Gunn-Pet erson 
trough in quasar absorption lines ( iGunn & PetersonI 1 19651) and 
in WMAP observat ions of the optical depth to recombination 
jSpergel et alj l2006h . Current observations of the temperature 
evolution of the IGM are similarly limited. At low redshift, 
observations of the Lya forest place constraints on the temperature 
of the I GM after reionization jSchave et a l. 2000: McDonald e t au 



20011: 



iZaldarriaga. Hui & TegmarkI 12001 : Theuns et al, 20021 : 

Hui&Haimad l2003h . Unfortunately, photoionization during 
reionization causes a large temperature increase that essentially 
erases information about the preceding period. At high redshift, it 
is assumed that the gas cools adiabatically after thermal decoupling 
from the C MB at z ~ 1 50, when Compton scattering becomes 
inefficient jPeeblesI [l993h . The intermediate regime, where the 



first sources have "switched on", is poorly constrained. Once 
collapsed structures form many diff erent heating mec h anism s 
are possible, e.g., shock heating jFurlemetto & Loea 2004 ). 
reson ant scattering of Lya ph otons ( Madau, Meiksin & Rees 
19971 : IChen & Miralda-Escud^ [2004 : Chuzhov & Shapirc 
20061 : iMeiksinI l2006l : iRvbickj l2006l: Furlanetto & Pritchard 

20011 



2006h. and X-ray heating jOstriker & Gn edin 1996 : lOh ik, . . 
Venkatesan. Giroux & Shull 1200 ll : iRicotti & Ostrikeil 120041) . 



Determining the thermal history and identifying the important 
heating mechanisms requires new observations. 

Future telescopes such as the James Webb Space Telescope 
hope to image high-redshift sources directly. However, seeing 
the sources is not the same as seeing the heating and ionization 
they cause in the IGM. The most promising technique for prob- 
ing the thermal history of the IGM before reionization is via ob- 
servation of the 21 cm hyperfine transition of neutral hydrogen 
jFurlanetto et alj|2006l and references therein). This line may be 
seen in absorption against the CMB, when the spin temperature 
Ts is less than the CMB temperature T.y, or in emission, when 
Ts > T^. Three prototype low-frequency interferometers (LO- 
FAfQ, MW^E and PASlO) are under construction and should be 
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capable of observing the redshifted 21 cm signal from gas at red- 
shifts z < 12, with the proposed Square Kilometre ArrajQ (SKA) 
capable of probing even higher redshifts. A great deal of theo- 
retical work has now been d one in calculating the 21 cm signal 



from fluctuations in density S i Loeb & Zaldarriaga"2004'), t he Lyg 
flux Jc, (Barkana & Loebl l2005bl: IPritchard & Furlanet tol l2006h . 
and the neutral fraction jZ aldarriaga, Furlanetto & Hernq ulsll2004l : 
iFurlanetto. Zaldarriaga & H ernquist 2004). Fluctuations in the 21 
cm brightness temperature Ti, also occur because of fluctuations in 
the gas kinetic temperature Tk, but this has not yet been explored. 

In this paper, we explore the effect of inhomogeneous X- 
ray heating by the first luminous sources on the 21 cm signal 
using analytic techniques. We first bu ild a model for th e global 
thermal history of the IGM following iFurlanettd j2006l) . In this 
model, we assume that a population of X-ray sources resulting 
from the remnants of the first s tars is responsible for heating the 
IGM dOstriker & Gnedin' ^igg^: [oi]|200ll : K^nkatesan et al.ll200ll : 
iRicotti & Ostriker ,2004.) . X-ray heating is dominated by soft X- 
rays (E < 2 keV), as harder X-rays have a mean free path com- 
parable with the Hubble scale. These long mean free paths have 
often motivated the simplifying assumption that X-rays heat the 
IGM uniformly. In fact, clustering of the sources and the 1 /r^ de- 
crease of flux with distance combine to produce significantly in- 
ho mogeneous heating. We develop a formalism, based upon that 
of iBarkana & Loea d2005b ). for calculating the temperature fluc- 
tuations that are sourced by these inhomogeneities. We use this to 
explore features in the 21 cm power spectrum that constrain the 
evolution of Tk ■ This calculation also motivates a consideration of 
the possibility of using 21 cm measurements to constrain the X-ray 
emission spectrum of the first sources. 

Simulations of the early universe have yet to address the spec- 
trum of temperature fluctuations in the period before reioniza- 
tion. Previous analytic consideration of fluctuations in Tk has fo- 
cused on_ttie_genodfoUowiM sources 
form tearkan a & Loebll2005 JiNaoz & Barkanall2005l) . Tempera- 
ture fluctuations induced by the first sources have not previously 
been considered in detail. 

The 21 cm signal can be thought of as a tool for prob- 
ing various radiation backgrounds. Gas temperature fluctuations 
probe the X-ray background, neutral fraction fluctuations probe 
the ionizing UV background, and Lya fluctuations probe the non- 
ionizing UV background. While the focus of this paper is X-ray 
heating of the IGM, in practice, the different sources of 21 cm 
fluctuation are not cleanly separated. In order to properly estab- 
lish context, we briefly re-examine the signal from fluctuations in 
the L yg flux, incorporating Lyg production by X-ray excitation 



^ Lyc 

of HI dChen & Miraida-Escudel2006l : [chuzhov, Alvarez & Shapirol 
l2006t) . and determine whether this contains extra useful informa- 
tion for constraining the spectral properties of the X-ray sources. 
Finally, we explore the feasibility of separating information on the 
temperature and Lyg flux fluctuations with the 21 cm signal. 

The layout of this paper is as follows. We begin by setting out 
the physics of the 21 cm signal in Sj2] Calculating this requires a 
model for the global history of the IGM, which we outline in ij3] 
Having established the mean history, in 21 we describe our frame- 
work for calculating fluctuations in Tk, Ja, and the neutral frac- 
tion. This is used to calculate the power spectrum for fluctuations 
in Tk in ij5] We then calculate the 21 cm signal in fjS] exploring 
the redshift evolution and dependence on the X-ray source spec- 



" See lCarilli & Rawlingj iiooj) . 



trum and luminosity. Finally, we discuss the possibility of obser- 
vationally detecting and separating these signals in iJ7]before con- 
cluding in ijs] Throughout this paper, we assume a cosmology with 
Qm = 0.26, = 0.74, Qt = 0.044, H = 100/ikms"^ Mpc"^ 
(with h = 0.74), ns = 0.95, and ag = 0.8, consistent with 
the latest measurements dSpergel et alj2006l) . although we have in- 
creased as above the best-fit WMAP value to improve agreement 
with weak-lensing data. 



2 21 CM SIGNAL 

We begin by briefly summarising the physics of the 21 cm signa l 
and refer the interested reader to IFurlanetto. Oh & Briggj d2006l) 
for further information. The 21 cm line of the hydrogen atom re- 
sults from hyperfine splitting of the IS ground state due to the in- 
teraction of the magnetic moments of the proton and the electron. 
The HI spin temperature Ts is defined via the number density of 
hydrogen atoms in the IS singlet and triplet levels, no and ni re- 
spectively, m/no = {gi/go)exp{-Ti,/Ts), where (ffi/ffo) = 3 
is the ratio of the spin degeneracy factors of the two levels, and 
T* = hc/k\2\cm = 0.0628 K. The optical depth of this transition 
is small at all relevant redshifts, so the brightness temperature of 
the CMB is 



n = 27xm{l + 5b) 



Q.bh'^ \ ( 0.15 1 + z\ ^''^ (Ts - T^ 



Qmz)\VL^h^ 10 y \ Ts ' ^^-^ 

Here xm is the neutral fraction of hydrogen and 5b is the fractional 
overdensity in baryons. The spin temperature is given by 



To 



(2) 



1 ~i~ Xct ~t~ Xc 

where Ta is the colour temperature of the Lyg radiation field at the 
Lyg frequency and is closely coupled to Tk by recoil during re- 
peated scattering. The spin temperature becomes strongly coupled 
to the gas temperature when xtot = Xc + Xa > 1. 
The collisional coupling coefficient is given by 

'^"^* [«:f-^(Tfc)nH + Ki^o(T'fc)ne] , (3) 

^ is the spontaneous emis- 
sion coefficient, K^-If) is tabulated as a function of Tk 



where Ai 



SAioTj 
= 2.85 X 10"^^ 



dAllis on & Dalgarno 1969; Zvgelman 2005) and kI_o is taken 
from Furlanetto & Furlanetto ( 2 0061). For a more d etailed analysis 

of the collisional coupling, see Hirata & SigurdsonI 1 2006). 

The Wouthysen-Field effect dWouthuvsenlll952|; |Field|[T95^ 
coupling is given by 



(4) 



27AioT-,meC 

where fa ~ 0.4162 is the oscillator strength of the Lyg transi- 
tion. Sa is a correction factor of order unity, which describes the 
detailed structure of t he photon distribution in the n ei ghbourhood 
of the Lyg resonanc e (Chen & Miralda-Escude 2004; Hirata 2006|; 
IChuzhov & Shapirdl2006 ; Furlanetto & Pritchard 2006). We make 
use of the approximation for Sa outlined in Furlanetto & Pritchard 
d2006h . For the models considered in this paper, Lyg coupling dom- 
inates over collisional coupling. 

Fluctuations in the 21 cm signal may be expanded 
dFurlanetto et alj2006l) 



5ti, = P5 + Px5x + I3a5a + I3t5t — 5;-)v 



(5) 
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where each 5i describes the fractional variation in the quantity i: 5a 
for fluctuations in the Lya coupling coefficient, 5x for the neutral 
fraction, St for Tk, and Sg^ for the line-of-sight peculiar velocity 
gradient. The expansion coefficients are given by 



/3 = 1 + 



/3x = 
Pa = 
Pt = 



1 + 



a:tot(l + xtot) ' 

(1 + ^tot ) 



(6) 



a;tot(l + Xtot) ' 

TK-Ty 
1 

r. T 



H d log kI" hh d log nfo^ 



dlogTK 



dlogTA- 



In this, we assume that baryons trace the density field exactly so 
that 5b = 5. All of these quantities are positive, with the excep- 
tion of Pt, whose sign is determined by {Tk — T-,). The appar- 
ent divergence in Pt when Tk = T7 is an artefact of expanding 
the fractional brightness temperature about a point where the mean 
brightness temperature fj, — 0. The physical quantity Ti,Pt is al- 
ways well behaved. 

Noting that in Fourier space Sg^ — —^^5 jBharadwai & Alii 

l2004h . where ^ is the angle between the line of sight and the 
wavevector k of the Fourier mode, we may us e equation ^ to form 
the power spectrum tearkana & Loebll2005ah 



Pt, (k, fi) = P^o (k) + (fc) + ^**P„4 (fc). 



(7) 



In theory, high precision measurements of the 3D power spectrum 
will allow the separation of these terms by their angular depen- 
dence. However, it is unclear whether the first generation of 21 
cm experiments will be a ble to achieve the high signal-to-noise re- 
quired for this separation iMcOuinn et al.l2005h . Instead, they will 
measure the angle averaged quantity 



Pt, (k) = P.a (k) + P^2 (fc)/3 + P„4 {k)/5. 



(8) 



In presenting our results, we will concentrate on P^j^2(k), which 
most cleanly separates out the different types of fluctuation, and 
PT,{k), which is easiest to observe. We will typically plot the 
power per logarithmic interval A = [k''^P(k) /2tt^]^^'^. 



3 GLOBAL HISTORY 
3.1 Outline 

We may express Tt as a function of four variables Tt = 
Ti){Tk ,Xi^ J a, uh)- In calculating the 21 cm signal, we require 
a model for the global evolution of and fluc tuations in t hese q uan- 
tities. We will follow the basic formalism of lFurlanettd ( l2006h , but 
first let us consider the main events in likely chronological order. 
This determines redshift intervals where the signal is dominated by 
fluctuations in the different quantities. 

z > 200: After recombination, Compton scattering maintains 
thermal coupling of the gas to the CMB, setting Tk ~ T^ so that 
we expect Tt — 0. 

40 < .z < 200: In this regime, adiabatic cooling means Tk < 
T~i and coUisional coupling sets Ts < T^,, leading to T), < 
and a possible absorption signal. At this time, Tt fluctuations are 
sourced by density fluctuations, potentially allowing cosrnology to 
be probed fcoeb & Zaldarriagall2004lHirata & Sigurdson"200g). 

^* ^ -2 ^1 40: As the expansion continues, decreasing the 



gas density, coUisional coupling becomes ineffective, absorption of 
CMB photons sets Ts = T-y, and there is no detectable 21 cm 
signal. 

z < Zi,: Once the first sources switch on at 2*, they 
emit both Lya photons and X-rays. In general, the emissivity re- 
quired for Lya coupling is significantly less than that for heating 
Tk above T^. Thus, in the simplest models, we expect the red- 
shift Za, where Lya coupling saturates Xa S> 1, to be greater 
than Zh, where Tk ~ T^. In this regime, T5 ~ < T-y and 
there is an absorption signal. Fluctuations are dominated by density 
fluctuations euid variation in t he Lya flux jSa rkana & Loeb 2005H; 
iPritchard & Furlanettdl2006l : IChen & Miraida-Escuddl2006l) . 

^ z < Za'- After Lya coupling saturates, fluctuations in 
the Lya flux no longer affect the 21 cm signal. By this point, heat- 
ing becomes significant and gas temperature fluctuations source Tt 
fluctuations. While Tk remains below T^ we see a 21 cm signal in 
absorption, but as Tk approaches T^ hotter regions may begin to 
be seen in emission. 

^ z < Zh'- After the heating transition, Tk > T-, and 
we expect to see a 21 cm signal in emission. The 21 cm bright- 
ness temperature is not yet saturated, which occurs at zt, when 
Ts ~ Tk 3> T-y. By this time, the ionization fraction has likely 
risen above the percent level. Brightness temperature fluctuations 
are sourced by a mixture of fluctuations in ionization, density and 
gas temperature. 

^ z < Zt'- Continued heating drives Tk 3> T-, at zt and 
temperature fluctuations become unimportant. Ts ~ Tk S> T7 
and the dependence on Ts may be neglected in equation (l^, 
which greatly simplifie s analysis of the 21 cm power spectrum 
(Santos & Coorav 2006). By this point, the filling fraction of HII 
regions probably becomes signific ant and ionization fluc tuations 
begin to dominate the 21 cm signal jFurlanetto et alj2004h . 

z < Zr'- After reionization, any remaining 21 cm signal origi- 
nates from overdense regions of collapsed neutral hydrogen. 

Most of these epochs are not sharply defined, so there should 
be considerable overlap between them. This seems the most likely 
sequence of events, althoug h there i s cons iderable uncertainty in 
the ordering of Za and Zh - iNussetl ( l2005h explores the possibil- 
ity that Zh > Za, so that X-ray preheating allows coUisional cou- 
pling to be important before the Lya fl ux becomes significant - 
Simulations of the very first mini-qua sar jKuhlen & Madaull2005l: 
iKuhlen. Madau & Montgomervll200^ also probe this regime and 
show that the first luminous X-ray sources can have a great im- 
pact on their surrounding environment. We note that these authors 
ignored Lya coupling, and that an X-ray background may gener- 
ate significant Lya photons jChen & Miralda-Escudell2006l) . as we 
discuss in i]3.3l 

In this paper, we will concentrate on the period after z*, when 
luminous sources "switch on", but before the IGM has been heated 
to temperatures Tk ^ T-j (our zt)- In this regime, Lya coupling 
dominates and the 21 cm signal is seen in absorption at high z but 
in emission at lower z. We shall explore this transition in more 
detail below. One of our key observables for 21 cm observations is 
the sign of Pt, which indicates whether Tk > T7 (provided that 
coUisional coupling can be neglected). 

3.2 Heating and ionization 

Having set the broad context, let us tighten our discussion with 
a concrete model for the evolution of the IGM; in this we follow 
iFurlanettol bOOeV We will distinguish between the ionization frac- 
tion Xi, relating to the volume filled by the highly ionized HII re- 
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gions that are located around clusters of sources, and the free elec- 
tron fraction Xe of the largely neutral gas outside these HII regions. 
The former is important for determining when reionization occurs, 
while the latter governs X-ray heating in the bulk of the IGM. We 
note that the volume filling fraction of the HII regions is well ap- 
proximated by Xi, which we will use to calculate volume averaged 
quantities. We further distinguish between Tk, the temperature of 
the IGM outside the HII regions, and the temperature of these pho- 
toionized regions Thii ~ 10* K. At high z, these regions are small 
and will not have a significant effect; while at low z, where reion- 
ization is well advanced, these HII regions will dominate and inval- 
idate our formalism. 

We begin by writing down equations for the evolution of Tk, 
Xi, and 

dT. 2r.dn 2 ^ 

.sA'D ^ — ^ n. 



3n At 
= (1 - a;e)Ai 
= {l-x^)A^ 



3fcs ^ — ' n 
j 

aACxluH , 



ctACx^riH , 



(10) 
(11) 



dt 

Axi 
'At 

AXe 

"dT 

where ej is the heating rate per unit volume, and we sum over all 
possible sources of heating/cooling j. We define Ai to be the rate 
of production of ionizing photons per unit time per baryon applied 
to HII regions, Ae is the equivalent quantity in the bulk of the IGM, 
ua = 4.2 X 10~^^cm^s~^ is the case-A recombination coeffi- 
cien{3at T = 10* K, and C = {nl) / {n^)'^ is the clumping factor. 
We model the clumping factor using C = 2; this value for C repro- 
duces the qualitative form of the histories inifurlanetto (2006) and 
ensures reionization at jz > 6. This approximation is appropriate 
only while Xi is small, and will fail towards the end of reioniza- 
tion, when clumpi ng becomes important in determining the effect 
of recombinations dMiralda-Escude. Haehnelt & Reesll2000h . 
In modelling the growth of HII regions, we take 



A,: = C(^) 



d/c, 



dt 



(12) 



where fcoii{z) is the fraction of gas inside collapsed objects at z 
and the ionization efficiency parameter C, is given by 



(13) 



with Mon the number of ionizing photons per baryon produced 
in stars, /* the fraction of baryons converted into stars, /osc the 
fraction of ionizing photons that escape the host halo, and Auc a 
correction factor for the presence of Helium. This model for Xi is 
motivat ed by a picture of HII r egions expanding into neutral hy- 
drogen iBarkana & Loebl [20011) . In calculating /coii, we use the 
IPress & Schechteil ( 1 19741) mass function An /Am and determine a 
minimum mass rrimin for collapse by requiring the virial temper- 
ature Tvir > 10* K, appropriate for cooling by atomic hydrogen. 
Decreasing this minimum mass, say to that of molecular cooling, 
will allow star formation to occur at earlier times shifting the fea- 
tures that we describe in redshift. We note that Xe ^ 1 at all red- 
shifts under consideration, as once the free electron fraction reaches 
a few percent further X-ray energy is deposited primarily as heat, 
not further ionization. 

To integrate equation we must spec ify which heating 
mechanisms are important. iFurlanettal j2006l) considers several 



^ Note that we use the case-A value, which amounts to assum- 
ing that ionizing pho tons are absorbed inside dense, neutral systems 
iMiralda-Escude et alJlioOOl) 



heatin g mechanisms including shock heating jpurlanetto & Loeb 



2004 ) and reso nant scattering of Lya photons ( Madau et aU 



1997: IChen & Miralda-Escude 2004; Chuzhov & Shapiro 2006 



Furlanetto & Pritchard 2006). We shall neglect these contributions 



to heating of the IGM, focusing instead on the dominant mecha- 
nisms of Compton heating and X-ray heating. While shock heating 
dominates the thermal balance at low z, during the epoch we are 
considering it, probably, heats the gas only slightly before X-ray 
heating dominates. 

Compton heating serves to couple Tk to T-y at redshifts z > 
150, but becomes ineffective below that redshift. In our context, it 
serves to set the initial conditions before star formation begins. The 
heating rate per particle for Compton heating is given by 



2 Ecompto 

3" 



^^{l + z)\ 



(14) 



1 + ,/"hc + Xe 1-7 u>^ 

where /hc is the helium fraction (by number), it^ is the energy 
density of the CMB, ar = 
cross-section, and we define 



t: 



.55 X 10" 



(15) 



3meC 

X-rays heat the gas primarily through photo-ionization of HI 
and Hel: this generates energetic photo-electrons, which dissipate 
their energy into heating, secondary ionizations, and atomic excita- 
tion. With this in mind, we calculate the total rate of energy depo- 
sition per unit volume as 



ex 



4?™,; / Au a,^^iJ,y{hu — hi/th), 



(16) 



where we sum over the species i —HI, Hel, and Hell, rii is the 
number density of species i, huth ~ Eth is the threshold energy 
for ionization, a^^i is the cross-section for photoionization, and J„ 
is the number flux of photons of frequency u. We may divide this 
energy into heating, ionization, and excitation by inserting the fac- 
tor fiiy), defined as the fraction of energy converted into form 
i at a specific frequency. The relevant division of the X-ray en- 
ergy depen ds on Xe and is c alculated using the fitting formulae of 
iShuU & van Steenber 3 11985.) . The fi{y) are approximately inde- 
pendent of V for hv > 100 eV, so that the ionization rate is related 
to the heating rate by a factor /ion/ (/hoat^^th). The X-ray number 
flux is found from 

Jx{z) = / Ai^Jxii',z), (17) 
nil 

where ix{v,z) is the comoving photon emissivity for X-ray 
sources, and i/' is the emission frequency at z' corresponding to 
an X-ray frequency at z 

., ..(1 + ^0 



(1 + -) 



(18) 



The optical depth is given by 



t{u,z,z')= f -^Az" [nmumiv") + nHcicTHciii^") 
Jz d^" 

+ nHoIlO-HcIl(!^")]i (19) 

where we calculate the cross-sections using the fits of IVemer et al.l 
( Il996h . Care must be taken here, as the cross-sections have a strong 
frequency dependence and the X-ray frequency can redshift consid- 
erably between emission and absorption. In practice, the abundance 
of Hell is negligible and may be neglected. 
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X-ray heating is often portrayed as uniform, as tlie X-ray pho- 
tons possess long mean free paths. The comoving mean free path 
of an X-ray with energy E is ( Furlanetto et al. 2006) 



A, 



-1/3 



15 



/ E 



V300 eV 



Mpc. (20) 



Thus, the Universe will be optically thick, over a Hubble length, 
to all photons with energy below E ^ 2[(1 + z)/15Y^'^xKi keV. 
The E~^ dependence of the cross-section means that heating is 
dominated by soft X-rays, which do fluctuate on small scales. In 
addition, though, there will be a uniform component to the heating 
from harder X-rays. 

We consider three possible sources of X-rays: starburst 
galaxies, supernova remnants (SNR), an d miniquasars ( Iohll200ll : 
iblover&B raiid' '20031: l^rlanitt^lSo^. The incidences of star- 
bursts and supe rnova remnants are like ly to be tied to the global star 
formation rate jGlover & Brandll2003h . For simplicity, we will as- 
sume that miniquasars similarly track the star formation rate (SFR). 

In reality, of course, their evolution could be considerably more 

I "1 1 

complex l IMadau et a l. 2004). We characterise these sources by an 

emissivity per unit (comoving) volume per unit frequency 



SFRD 



(21) 



where SFRD is the star formation rate density, and the spectral dis 
tribution function is a power law with index as 



(22) 



and the pivot energy hi^o — 1 keV. We assume emission within the 
band 0.1-30 keV, and set Lo = 3.4 x 10*° /x ergs"^ Mpc~^ 
where fx is a highly uncertain constant factor. This normalisation 
is chosen so that, with fx~l, the total X-ray luminosity per unit 
SFR is consistent with that observ ed in starburst galaxies in the 
present epoch (see lFurlanettdl20o3 for further details). Extrapolat- 
ing observations from the present day to high redshift is fraught 
with uncertainty, and we note that this normalisation is very uncer- 
tain. The total X-ray luminosity at high redshift is constrained by 
observations of the present day soft X-ray back ground, which rules 
out complete rei onization by X-rays, but allows considerable lati- 
tude for heating toiikstra, Haiman & Loebll200 4V Similarly, there 
is significant uncertainty in the spectra of these objects. We choose 
as = 1.5 f or starbursts, as = 1-0 for SNR, and as = 0.5 for 
miniquasars jMadau et al]|2004l) . These span the reasonable spec- 
tral dependence of possible X-ray sources. 

As in equation ilH . we model the star formation rate as track- 
ing the collapse of matter, so that we may write the star formation 
rate per (comoving) unit volume 



SFRD = p°(z)f*-^^fcou{z). 



(23) 



where p|; is the cosmic mean baryon density today. This formalism 
is appropriate for z > 10, as at later times star formation as a result 
of mergers becomes important. 

3.3 Lya flux 

Finally, we must describe the evolution of the Lya flux. This is 
produced by stellar emission (Ja,*) and by X-ray excitation of 
HI (Ja,x). Photons emitted by stars, between Lya and the Ly- 
man limit, will redshift until they enter a Lyman series resonance. 
Subsequently, they may generate Lya photons, as discussed in 



IPritchard & Furlanetto! hOOdt) and Iffiratal ilOOdl) . The Lya flux 
from stars J^,* arises from a sum over the Lyn levels, with the 
maximum n that contributes rimax ~ 23 determi ned by the size of 
the HI I region of a typical (isolated) galaxy (see lBarkana & Loea 
for details). The average Lya background is then 



^ ^ /"recycle (^) 



(24) 



, {1 + zf 

dz 



47r H(z') 



where Zmax(w) is the maximum redshift from which emitted pho- 
tons will redshift into the level n Lyman resonance, v'^ is the emis- 
sion frequency at z corresponding to absorption by the level n at z, 
/recycle (ti) is the probability of producing a Lya photon by cascade 
from level n, and e, {y, z) is the comoving photon emissivity for 
stellar sources. We calculate e* (y, z) in the same way as for X-rays 
(eq. 21), and define e*(j^) to be the spectral distribution function of 
the stellar sources. We consider models with Pop. II and very mas- 
sive Pop. Ill stars. In each case, we take (u) to be a broken power 
law with one index describing emission between Lya and Ly/3, and 
a se cond describing emissi on between Ly/3 and the Lyman limit 
(see lPritchard & Furlanettj 2006 for details). 

Photoionization of HI or Hel by X-rays may also lead 
to the production of Lya photons. In this case, some of the 
primary photo-electron's ene rgy ends up in excitations of HI 
dShull & van Ste enberg* 'l98 ^, which on relaxation may generate 
Lyg photons I Madau et al. ll997l : IChen & Miralda-Escud^ l2006l : 
IChuzhov et al.l 20061) . This Lya flux Ja,x may be related to the 



X-ray heating rate as follows. The rate at which X-ray energy is 
converted into Lya photons is given by 



ex,, 



ex, heat 



/e: 



/"heat 



(25) 



where /ex and /heat are the fraction of X-ray energy going into ex- 
citation and heating respectively, and pa is the fraction of excitation 
energy that goes into Lya photons. We then find the Lya flux by 
assuming that this injection rate is balanced by photons redshifting 
out of the Lya resonance, so 

1 



c ex,, 



(26) 



47r hfa HVa 

IShull & van Steenber i ( Il985l) calculated /ex and /heat, but 
their Monte Carlo method, gives only a little insight into the value 
of Pa. Although excitations to the 2P level will always generate 
Lya photons, only some fraction of excitations to other levels will 
lead to Lya generating casc ades. The rest will end with two pho- 
ton decay from the 2S level. IShuU & van Steenbergl ( ll985h consid- 
ered a simplified atomic system, in which collisional excitations to 
n > 3 levels were incorporated by multiplyi ng the excita tion cross- 
section to the n = 2 level by a factor of 1.35 (IShuIlll979h . Thus, we 
might expect of order 0.35/1.35 ~ 0.26 of collisional excitations 
to end at an n > 3 level. 

We estimate pa by calculating the probability that a sec- 
ondary electron of energy E^cc will excite HI from the ground 
sta te to the l evel nL, using the collisional cross-section|f| 
of iBrav et alj ( |2002|) . and then applying the probability that 



Taken from |http://atom.murdoch.edu.au/CCC-WWW/index.htnJ] 
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Figure 1. Mean IGM thermal history for model A (thick curves) and B 
(thin curves), (a): Tk (solid curve), T-y (dashed curve), and Tg (dotted 
curve), (h): Volume averaged fj, (sohd curve). The zero line is indicated by 
a dashed horizontal line. Note that this is the thermal history outside of the 
ionized HII regions. 



Figure 2. Ionization histories for model A (thick curves) and B (thin 
curves), (a): Xi (dotted curve), Xe (dashed curve), and the volume averaged 
ionization fraction Xi = Xi + {I — Xi)xe (solid curve), (b): The quantities 
ftTf,. We plot 13 (solid curve), /Bx (long dashed curve, overlapping with /3), 
f3x (short dashed curve), and 13a (dotted curve). 



the resulting cascade will p rod uce a Lvq p hoton, taken from 
IPritchard & Furlanettd (120061) andlffiratal OOOeh . The iterative pro- 



cedure of IPritchard & Furlanettd feooeh gives the probability of 
producing a Lya photon by cascade from the level nL as: (0, 1) 
for (2S, 2P), (1, 0, 1) for (3S, 3P, 3D), and (0.584, 0.261, 0.746, 1) 
for (4S, 4P 4D, 4F). 

Summing over atomic levels n < 4, we obtain pa = 0.79 for 
Escc = 30 eV. The contribution from n > 4 levels is small as the 
coUisional cross-sections drop off rapidly as n increases. The exact 
result depends upon the energy distribution of the secondary elec- 
trons, which in turn depends upon the spectrum of ionizing X-rays. 
Our chosen value for Sscc corresp onds to the m ean electron energy 
(obtained using the distribution o f IShullI 19791) produced by X-rays 
of energy 1.7 keV, which is the mean X-ray energy from a source 
with spectral index a = 1.5 over the band 0.1 — 30keV. Calcu- 
lating Pa exactly requires an update of the IShull & van Steenberd 
i ll 985b calculation, but, by considering different values for Escc, we 
conclude that it should differ frompa ~ 0.79 by less than 10%. 



3.4 Model histories 

Having outlined the various elements of our global history, we will 
restrict ourselves to considering two models. These will be A) Pop. 
II stars -I- starburst galaxies and B) Pop. Ill -I- starburst galaxies. 
Of course, these are only two of an infinite set of possibilities, but 
they serve to illustrate the effect of different Lya and X-ray lumi- 
nosities on the signal. We use parameters corresponding to Pop. II 
(/esc = 0.1, /* = 0.1, iVion = 4000) and ver y massive Pop. Il l 
(Uc = 0.1, = 0.01, iVion = 30000) stars ( lFurlanettdl20o3) . 
although we note that these values are highly uncertain. We take 
/jf = 1 in both models, to allow straightforward comparison be- 
tween the two models. The amplitude of the X-ray background is 



extremely uncertain, so that fx is almost unconstrained, and we 
defer discussion of its effects until %.3\ 

The global histories produced by these models are shown in 
Figure[T]and Figure|2] In Figure[T] we see the evolution of Tk, Ts, 
and T-y . Note that, while both models produce the same qualitative 
behaviour, the reduced star formation rate in Model B delays the 
onset of heating from z ~ 18 to z ~ 15. We also see that the heat- 
ing transition, where Tk ~ T7, occurs at Zh ~ 14 in Model A and 
Zh ~ 11 in Model B. We have assumed that the X-ray luminos- 
ity per unit star formation is the same for both populations, so this 
is an effect of /* rather than the initial mass function. In addition. 
Pop. Ill stars produce fewer Lya photons than Pop. II stars further 
slowing the onset of Lyman coupling. 

Figure [ij) shows the distinctive Tj, signature of absorption at 
z > Zh followed by emission at z < Zh in both models. The 
signal is signific antly larger and more extended in Model A (See 
lFurlanettdl2003 for more detailed discussion of such histories). 
The ionization history is outlined in Figure [2^ and shows that Xi 
evolves similarly in both models, as they have similar values for 
^. The electron fraction in the IGM Xe is depressed in model B, 
where there is a smaller X-ray background. Note that Xc remains 
much smaller than Xi once ionization begins. Both ionization his- 
tories produce an optical depth to the surface of last scattering 
Tri ~ 0.07, consistent with the WMAP t hird year observations 
of Tri = 0.09 ± 0.03 dSpergel et al]|2006i) . although slightly on 
the low side. Our model for temperature fluctuations will be geared 
towards making predictions for the largely neutral IGM outside of 
the ionized HII regions surrounding clusters of UV sources. Con- 
sequently, from Figure |2^, we expect our model to be valid for 
z > 12, where Xe < 0.1 and the filling fraction of the HII regions 
is small. 

Figure|2j5 shows l3iTt, which is a measure of the sensitivity of 
the 21 cm signal to fluctuations in each fundamental quantity. If the 
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21 cm signal were dominated by component i and if the fluctuation 
had unit amplitude 5^ ~ 1, then PiTb gives the amplitude of the 
21 cm signal. Note that the curves for 5 and 5x are almost indistin- 
guishable and track Tj. In contrast, the curves for Sa and 5t show 
clear peaks - representing windows where an existing signal might 
be seen. We may identify Zh as the point where Ti, = and all 
curves except that for /3t go to zero. At this point, the only fluctua- 
tions in Tb arise from fluctuations in Tk- In practice, this "null" is 
more mathematical than physical, as inhomogeneities will blur the 
situation. The redshift window for observing the 21 cm signal is 
clearly much narrower in model B, indicating that it will be much 
more confused than in model A. 



4 FORMALISM FOR TEMPERATURE AND 
IONIZATION FLUCTUATIONS 

Having specified our global history, we now turn to calculating the 
fractional fluctuations 5a, St, and S^- Note that we will primar- 
ily be interested in the signal from the bulk of the IGM, working 
at redshifts where Xi < 0.1, so that we will ignore the fluctua- 
tions induced by HII regions. We begin by forming equations for 
the evolution of St and Se (the fractional fluctuation in a: e) by per- 
turbing equations lO and dl lb (see also ISarkana & Loebll20053 : 
iNaoz & Barkanall2005l) . This gives 



dST 
_ (1- 



2d5_ 



E 



2Ahe 



'-[Sa 



(27) 



-A,[5Ae - 5e] - aACx^nnlSe + S], (28) 



dSe 

dt 

where an overbar denotes the mean value of that quantity, and 
A = e/n is the ionization or heating rate per baryon. We also need 
the fluctuation in the neutral fraction Sx = —x^/il — Xs)Se and in 
the Lya coupling coefficie nt Sg = Sj^ , neglecting the m ild tem- 
perature dependence of Sa dFurlanetto & Pritchardl l20061. 

To obtain a closed set of equations, we must calculate the fluc- 
tuations in the heating and ionizing rates. Perturbing equation l ll4t 
we find that the contribution of C ompton scattering to th e right 
hand side of equation Ml\ becomes jNaoz & BarkanalboOSi) 



2Ah 



1 + /ho + a; 



Tk 



1 S' 



St) 



(29) 



where St-^ is the fractional fluctuation in the CMB temperature, 
and we have ignored the effect of ionization variations in the neu- 
tral fraction outside of the ionized bubbles, which are small. Before 
recombination, tight coupling sets Tk = T7 and St = St^ . This 
coupling leaves a scale dependent imprint in the temperature fluc- 
tuations, which slowly decreases in time. We will ignore this effect, 
as it is small (~ 10%) below 2 = 20 and once X-ray heating be- 
comes effective any memory of these early temperature fluctuations 
is erased. At low z, the amplitude of <5t-, becomes negligible, and 
equation i29\ simplifies. 

Our main challenge then is to calculate the fluctuations in the 
X -ray heating. We s hall a chieve this by paralleling the approach 
of iBarkana & Loebl JioOSbi) to calculating fluctuations in the Lya 
flux from a p opulation of stellar sources. We first outline their re- 
sults (see also lPritchard & Furlanettoll2006 ). Density perturbations 
at redshift z' source fluctuations in Ja seen by a gas element at 



redshift z via three effects. First, the number of galaxies traces, 
but is biased with respect to, the underlying density field. As a re- 
sult an overdense region will contain a factor [1 + b{z')S] more 
sources, where b{z') is the (mass-averaged) bias, and will emit 
more strongly. Next, photon trajectories near an overdense region 
are modified by gravitational lensing, increasing the effective area 
by a factor (1 + 25/3). Finally, peculiar velocities associated with 
gas flowing into overdense regions establish an anisotropic redshift 
distortion, which modifies the width of the region contributing to a 
given observed frequency. Given these three effects, we can write 
Sa = Sj^ — Wa{k)S, where we compute the window function 
M^Q,*(fc) for a gas element at z by adding the coupling due to Lya 
flux fr om each of the Lyrt res onances and integrating over radial 
shells l lBarkana & Loebl(2005bl) 



WaAk) = 



-E 

a,* „ 

71 — 2 

D{z) 



az — ; 

Az' 



[l + fe(2')bo(fcr)-|j2(fcr)|, (30) 



where D{z) is the linear growth function, r — r{z, z') is the dis- 
tance to the source, and the ji (x) are spherical Bessel functions of 
order I. The first term in brackets accounts for galaxy bias while the 
second describes velocity effects. The ratio D{z') / D{z) accounts 
for the growth of perturbations between z' and z. Each resonance 
contributes a differential comoving Lya flux d ji"' / dz', calculated 
from equation i24l . 

We plot Wa,*{k) in Figure[3] On large scales, Wa,*{k) ap- 
proaches the average bias of sources, while on small scales it dies 
away rapidly encoding the property that the Lya flux becomes 
more uniform. In addition to the fluctuations in Ja.*, there will 
be fluctuations in Ja.x- We calculate these below, but note in 
passing that the effective value of Wa is the weighted average 
Wa ~ Wa,i{Ja,i/Ja) of the Contribution from stars and X- 
rays. 

We now extend the formalism of iBarkana & Loebl ( l2005bh 
in an obvious way to calculate fluctuations in the X-ray heating 
rate. First, note that for X-rays (Jaj^,, = <5Ai,cat ~ ~ ^^x ' 
as the rate of heating, ionization, and production of Lya photons 
differ only by constant multiplicative factors (provided that we 
may neglect fluctuations in Xe, which are small). In each case, 
fluctuations arise from variation in the X-ray flux. We then write 
Sax = Wx{k)S and obtain 



Wx 



Ax JEtu Jz 



j dkxjE) 
' dz' 



Djz') 
D{z) 



[l + &(2')bo(fcr)-|j2(fcr)|, (31) 



where the contribution to the energy deposition rate by X-rays of 
energy E emitted with energy E' from between redshifts z' and 
z' + dz' is given by 



dAxjE) 47V ,^^ dJx{E,z) ^^ 

^ T'^-^^) dz' - 



(32) 



and Ax is obtained by performing the energy and redshift inte- 
grals. Note that rather than having a sum over discrete levels, as 
in the Lya case, we must integrate over the X-ray energies. The 
differential X-ray number flux is found from equation 

The window function Wx (k) gives us a "mask" to relate fluc- 
tuations to the density field; its scale dependence means that it is 
more than a simple bias. The typical sphere of influence of the 
sources extends to several Mpc. On scales smaller than this, the 
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Figure 3. Wa,i,(k) (dotted curves), Wx{k) (dashed curves), and gT(k) 
(solid curves) at 2 = 20 (thin curves) and z = 15 (thick curves) for Model 
A. 



shape of Wx(k) will be determined by the details of the X-ray 
source spectrum and the heating cross-section. On larger scales, 
the details of the heated regions remain unresolved so that Wx{li) 
will trace the density fluctuations. 

A further word of explanation about this calculation is worth- 
while. An X-ray is emitted with energy E' at a redshift z' and red- 
shifts to an energy E at redshift z, where it is absorbed. To calculate 
Wx we perform two integrals in order to capture the contribution 
of all X-rays produced by sources at redshifts z' > z. The integral 
over z' counts X-rays emitted at all redshifts z' > z which redshift 
to an energy Eatz; the integral over E then accounts for all the X- 
rays of different energies arriving at the gas element. Together these 
integrals account for the full X-ray emission history and source 
distribution. Many of these X-rays have travelled considerable dis- 
tances before being absorbed. The effect of the intervening gas is 
accounted for by the optical depth term in Jx- Soft X-rays have 
a short mean free path and so are absorbed close to the source; 
hard X-rays will travel further, redshifting as they go, before being 
absorbed. Correctly accounting for this redshifting when calculat- 
ing the optical depth is vital as the absorption cross-section shows 
strong frequency dependence. In our model, heating is dominated 
by soft X-rays, from nearby sources, although the contribution of 
harder X-rays from more distant sources can not be neglected. 

We compare the form of Wx (k) and the stellar component 
of Wa{k) in Figure [3] Including the X-ray contribution in Wa{k) 
drives that curve towards the Wx (k) curve. Note that Wx shows 
significantly more power on smaller scales than Wa , reflecting the 
greater non-uniformity in the X-ray heating; most heating comes 
from soft X-rays, which have mean free paths much smaller than 
the effective horizon of Lya photons. Also, while Wa shows a sub- 
tle break in slope at A; ~ 3 Mpc^^, Wx shows no obvious features 
indicative of preferred scales. Both Wx and Wa trace the bias on 
very large scales. 

Returning now to the calculation of temperature fluctuations, 
to obtain solutions for equations \2H and ([28}, we let St = 



gT{k,z)5, Se = ge{k,z)S, 5a = Wa(k,z)5, and S a^ = 

Wx(k, z)S, foll owing the approach of Bharad wai & Alii ( |2004) . 
Unlike iBharadw ai & Ah (2004), we do not assume these quanti- 
ties to be independent of scale, and so we must solve the resulting 
equations for each value of k. Note that we do no t include the scale 
dependence induced by coupling to the CMB l lNaoz & Barkanal 
In the matter dominated limit, we have S cc (1 + z) ^ and 
so obtain 



gr - 2/3 



dz 

dz \1 + z 
where we define 



-Qx(z)[Wx{k)-gT]-Qc{z)gT, (33) 
Qi (z) [Wx (fc) - gel + Qr iz)[l + ge] , (34) 

(1 - Xe) Aion,X 



_ aACXeUH 



{l + z)H{z)' 



Qciz) 



(1 + z f 



1 + /ho + Xe U,H{Z) Tk ' 



and 



2Ah 



ZkBTK{l + z)H{z) 



(35) 



(36) 



(37) 



(38) 



These are defined so that Qr and Qi give the fractional change 
in Xe per Hubble time as a result of recombination and ionization 
respectively. Similarly, Qc and Qx give the fractional change in 
Tk per Hubble time as a result of Compton and X-ray heating. 
Immediately after recombination Qc is large, but it becomes neg- 
ligible once Compton heating becomes ineffective at 2 ~ 150. The 
Qr term becomes important only towards the end of reionization, 
when recombinations in clumpy regions slows the expansion of HII 
regions. Only the Qx and Qi terms are relevant immediately after 
sources switch on. We must integrate these equations to calculate 
the temperature and ionization fluctuations at a given redshift and 
for a given value of k. 

These equations illuminate the effect of heating. First, con- 
sider gr, which we can easily relate to the adiabatic index of the 
gas 7a by (7t = 7a — 1, giving it a simple physical interpretation. 
Adiabatic expansion and cooling tends to drive gr 2/3 (corre- 
sponding to 7a = 5/3, appropriate for a monoatomic ideal gas), but 
when Compton heating is effective at high z, it deposits an equal 
amount of heat per particle, driving the gas towards isothermality 
(gr 0). At low z, where X-ray heating of the gas becomes sig- 
nificant, the temperature fluctuations are dominated by spatial vari- 
ation in the heating rate (gr Wx)- This embodies the higher 
temperatures closer to clustered sources of X-ray emission. If the 
heating rate is uniform Wx{k) ~ 0, then the spatially constant 
input of energy drives the gas towards isothermality gr 0. 

The behaviour of g^ is similarly straightforward to interpret. 
At high redshift, when the IGM is dense and largely neutral, the 
ionization fraction is dominated by the recombination rate, pushing 
gx —1, because denser regions recombine more quickly. As the 
density decreases and recombination becomes ineffective, the first 
term of equation i34\ slowly pushes g^ — > 0. Again, once ioniza- 
tion becomes important, the ionization fraction is pushed towards 
tracking spatial variation in the ionization rate (gx Wx)- Note 
that, because the ionization fraction in the bulk remains less than 
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k [Mpc-'] 

Figure 4. Fluctuations in T[f. In each panel we plot TfcAj'(k) for the case 
of inhomogeneous X-ray heating (thick curves) at 2: = 20 (long dashed 
curve), 2 = 15 (dotted curve), 2: = 13 (short dashed curve), and 2 = 10 
(solid curve). For comparison, we plot the case of uniform heating al z = 
10 (thin solid curve) and 2 = 20 (thin long dashed curve), (a): Model A. 
(b): Model B. 



a few percent, fluctuations in the neutral fraction remain negligibly 
small at all times. 

The scale dependence of gx is illustrated in Figure [3] gr tries 
to track the heating fluctuations Wx{k) (as in the 2 = 15 curve), 
but two factors prevent this. First, until heating is significant, the 
effect of adiabatic expansion tends to smooth out variations in gr- 
Second, gr responds to the integrated history of the heating fluctu- 
ations, so that it tends to lag Wx somewhat. When the bulk of star 
formation has occurred recently, as when the star formation rate is 
increasing with time, then there is little lag between gr and Wx ■ 
In contrast, when the star formation rate has reached a plateau or 
is decreasing the bulk of the X-ray flux originates from noticeably 
higher z and so gr tends to track the value of Wx at this higher 
redshift. On small scales, the heating fluctuations are negligible and 
gr returns to the value of the (scale independent) uniform heating 
case. 



5 TEMPERATURE FLUCTUATIONS 

Before calculating the 21 cm signal, let us first examine the gas 
temperature fluctuations themselves. Figure |4] shows the power 
spectrum of temperature fluctuations Pt (fc) for models A and B re- 
spectively. We see that in both cases the fluctuations are small until 
z < 20. At lower redshifts and on larger scales (fc ~ 0.1 Mpc~^), 
the heating fluctuations source a significant (factor of ~ 50) en- 
hancement over the uniform heating case. This is to be expected. 
Uniform heating of the gas tends to erase temperature fluctuations, 
while inhomogeneous heating causes them to grow. Thus we ob- 
serve a huge increase in power. The fluctuation amplitude in Model 
B is generally smaller than in Model A as a consequence of the re- 
duced heating from the decreased SFR in Model B. In both cases, 
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Figure 5. Evolution of the fluctuations in Tx with redshift for Model A 
(thick curves) and B (thin curves), (a): We plot the amplitude of fluctuations 
in given by Ay (fc) at fc = 0.1 Mpc^ ^ in the case of uniform heat- 
ing (dashed curves) and when fluctuations in the heating rate are considered 
(solid curves). For comparison, we plot Ifj. — T-y | (dotted curves). Only in a 
small region of width Az 1 around Tk = Try do the fluctuations exceed 
this threshold, (h): Evolution of g^. We plot at fc = 0.1 Mpc~^ for 
the uniform (short dashed curves) and fluctuating cases (solid curves). We 
also plot Wx (dotted curve) and, for comparison, Wa^i,{k) (long dashed 
curve). Notice how qt rises to track Wx once heating becomes effective. 



the temperature fluctuations remain small, St < 1 (compare with 
Figure[Tll, justifying our linear approximations. 

Figure [5] illustrates the redshift evolution of the tempera- 
ture fluctuations. We choose to follow a single wavenumber k — 
0.1 Mpc~^, which is both within those scales accessible to future 
experiments and demonstrative of the effect. If the gas is heated 
uniformly (dashed curves), then gr rapidly becomes negligible 
once heating becomes effective. By depositing the same amount of 
energy per particle the gas is driven towards isothermality. When 
heating fluctuations are taken into account gr may grow or de- 
crease depending on scale. We observe that, for the scale chosen 
here, the amplitude of the temperature fluctuations grows steadily 
with time, but gr decreases. This is a consequence of the sources 
becoming less biased with time so that Wx{z) decreases with z. 
On very small scales, where Wx (k) is negligible, gr will trace the 
uniform heating curve. 

Recall that whether we observe the 21 cm line in emission 
or absorption depends on the sign of Tg ~ T-,. Assuming that 
Ts ~ Tk, when Tk < T^-, hotter regions have a spin temperature 
closer to the CMB temperature and so appear more faintly in ab- 
sorption. As heating continues, it is these regions that are first seen 
in emission, when their temperature exceeds T-y. Once Tk > T-,, 
these hotter regions produce the largest emission signal. 

We see from Figure [5] that for a short window around Zh 
(where Tk ~ T-,) temperature fluctuations may raise Tk above 
T-y in these hot regions, even when Tk is less than T^,. We inter- 
pret this to mean that within this window the 21 cm signal will be 
a confusing mix of emission, from hotter regions, and absorption, 
from cooler regions. In the case of uniform heating this window 
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Figure 6. Evolution of brightness temperature fluctuations for Model A 
(thick curves) and B (thin curves), (a): We plot |Ti,| Atj, (fc) at fc = 
0.1 Mpc"'^ including the effects of heating fluctuations (dotted curves), 
Ly« fluctuations (dashed curves), and both heating and Lya fluctuations 
(solid curves), (b): We plot If^p A^j {^) with the same line conventions. 



Figure 7. Full T;, power spectra for Model A. We plot the power spectra at 
2 = 20 (long dashed curve), z = 19 (short dashed curve), z = 18 (dotted 
curve), and z = 17 (solid curve), (a): |f(, | Atj_ (fc). We plot |fi,|Ai5 (thin 
solid curve) at z = 19 for comparison, (h): |Ti,p A^j C^^)- The sign of 
A^2 C^) is indicated as positive (thick curves) or negative (thin curves). 



is very narrow, but when fluctuations are included it extends to a 
significant (Az ~ 1) width. This indicates that the transition from 
absorption to emission will not be abrupt, but extended. 



6 21 CM POWER SPECTRUM 
6.1 Redshift evolution 

Finally, we write the full 21 cm power spectrum as 

PT,{k,^l)=fi{|3' +^^fPss{k), (39) 

where 

P' = 13- I3^x,g,/{1 + X,) + Prgr + Pc^W^- (40) 

Within our model we may neglect the term corresponding to the 
neutral fraction, as the free electron fraction in the IGM remains 
small at all times. We now consider how the 21 cm power spectrum 
evolves with redshift. 

Figure [6] shows the evolution of the brightness temperature 
fluctuations at a single scale k = 0.1 Mpc~^ with redshift. First, 
note that in the bottom panel A^2 changes sign when we include 
temperature fluctuations (note that A^2 is not an auto-correlation 
and so is free to have a negative sign). Physically, this occurs be- 
cause when Tk < there is an anti-correlation between Tt and 
Tk, i.e. increasing Tk decreases Tt- Observing P^2 < is a clear 
sign that Tk < T^- Mathematically, this can be seen because Pt is 
the only one of the fluctuation coefficients that can become nega- 
tive. Of course, if Pst or other cross-correlations become negative 
we can also get P^2 < 0, but this should not be the case for ra- 
diative heating or Lya coupling, as we expect emitting sources to 
be most common in overdense regions. Only in the case of P^s 
might we expect a negative cross-correlation, as increasing the UV 



radiation is likely to decrease the neutral fraction. In the high red- 
shift regime, before significant ionization has occurred, this term is 
negligible. 

Adding the Lya fluctuations, we see a clear double peaked 
temporal structure in the evolution of At^ , which is dominated by 
Lya fluctuations at high z and temperature fluctuations at lower 
z (were we to include the effects of ionization fluctuations, there 
would be a third peak at still lower redshift). We note that there is 
considerable overlap between the two signals, which will compli- 
cate extracting astrophysical information. The situation is similar 
in Model B, although here the relevant signal is compressed into 
a narrower redshift window. We note that the amplitude of fluctu- 
ations induced by the gas temperature is significantly larger than 
those from the Lya signal and present at lower redshifts. Both of 
these features make the temperature fluctuation signal a plausible 
target for future observations. 

To illustrate the scale dependence of this signal, we exam- 
ine a series of redshift slices. We will make plots for model A. 
Although the same evolution applies for model B, the events are 
shifted to lower redshift Az ~ 3 and the transitions are somewhat 
compressed in redshift. We begin by examining the high redshift 
regime, where Lya fluctuations dominate the 21 cm signal, but 
temperature fluctuations become important as we move to lower 
redshift. 

Figure |7] shows redshift slices from z = 17 — 20. We can 
see from Figure |6] that Lya fluctuations dominate the signal for 
z > 18. The z = 20 and z = 19 curves show the expected excess 
of p ower on large scales for Lya fluctuations from the first sources 
(see iBarkana & Loebll2005b for a full analysis of this signal). At 
z = 18, we begin to see the effects of the temperature fluctuations 
through the dip in power between fc = 0.1 and 1 Mpc~^. This dip 
occurs because Pt < 0, contrasting with the other Pi, which are 
positive. Physically, in this regime Tk < T^ and regions that are 
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Figure 8. Full Tj, power spectra for Model A. We plot the power spectra 
at 2 = 16 (solid curve), z = 15 (short dashed curve), 2 = 14 (dotted 
curve), and 2 = 13 (long dashed curve). Note that the 2 = 13 curve 
would be seen in emission, and the 2 = 14 curve in a mixture of emission 
and absorption. The other curves would be seen in absorption against the 
CMB. (a): \%\AT^(k). We plot \ (thin solid curve) at z = 14 for 

compaiison. (b): |f(,p {^)- The sign of Q') is indicated as positive 
(thick curves) or negative (thin curves). 



hotter have a smaller brightness temperature. In our model, denser 
regions are more strongly coupled, which increases Ti,, but are also 
hotter, which tends to decrease Tb- These two effects compete with 
one another and produce the dip. 

At 2: = 17, temperature fluctuations grow large enough to 
drive P' negative over a range of scales, where they outweigh the 
Lya fluctuations. This leads to a sign change in A^2, but also 
imprints a distinctive trough-peak-trough structure in Ay^. Here 
Lyo! fluctuations dominate on the largest scales, temperature fluc- 
tuations on intermediate scales, and density fluctuations on small 
scales. For this to occur, we require that > gr on large scales, 
which can only occur if Wa and gr show different scale depen- 
dence. This always occurs at some redshift in our model, as both 
Wa and Wx tend towards the same value on large scales, but gr 
lags behind (and so is smaller than) Wx on those scales. 

From Figure |6l we see that Tk fluctuations dominate at z < 
17 and that Lya fluctuations become negligible for 2 < 15. In Fig- 
ure[8]we plot redshift slices in the range z = 13 — 16. At 2 = 16 
and z = 15, we see a sign change in A^2 , which is a distinctive sig- 
nature of the temperature fluctuations when Tk < T-y . This is seen 
in Atj as a peak on large scales, followed by a trough at smaller 
scales. The position of the peak depends upon the shape of gr and 
thus the X-ray source spectrum. We will consider this in more detail 
in the next section. 

Notice that the heating transition occurs very close to 2 = 14, 
so that the 21 cm signal at this redshift would likely be seen in a 
mixture of absorption and emission. In addition, this curve is domi- 
nated by gas temperature fluctuations. We see this in Figure[8]where 
the contribution from density fluctuations at 2 = 14 (thin solid 
curve) is at least a factor of two smaller than At,, on all scales. Re- 



Figure 9. Effect of X-ray spectra on 2 1 cm power spectra. We show results 
at 2 = 15 for model A and take a = 1.5 (dotted curve), a = 1.0 (solid 
curve), and a = 0.5 (dashed curve), (a): \Ti,\Aj'^{k). We illustrate the 
uniform heating case by the thin solid curve, (b): |T[,pA^2 C^)- 



call from Figure|2]that when Ti, ~ only the combination |T(,|/3t 
is significant. 

The position of the sign change moves to smaller scales as 
the gas is heated and the temperature fluctuations become larger. 
Eventually, the IGM heats to Tk > T^, hotter regions have a higher 
brightness temperature than average and /3t > 0. Once this occurs 
the trough disappears entirely and the peak on large scales is no 
longer quite so distinctive (see 2 = 14 curve). The continued IGM 
heating drives /3t ^ and diminishes the effect of the temperature 
fluctuations. By 2 = 13 there is no longer a clear peak in either 
A^2 or Atj, , although there is still considerable excess power on 
large scales. By 2 = 10, Tk 3> T-, and temperature fluctuations 
no longer impact the 21 cm signal significantly. 

Once the ionization fraction becomes large (xi > 0.1), the 
21 cm signal becomes d ominated by the impr int of HII regions 
jZaldarriaga et al.l |2004| ; iFurlanetto et alj |2004l). This eventually 
produces a distinct knee in the 21 cm power spectrum resulting 
from the characteristic size of the bubbles. We note that our models 
have Xi < 0.1 at 2 > 12, so that we do not expect ionization fluctu- 
ations to significantly affect the results we have outlined for Model 
A. In the case of model B, temperature fluctuations remain signif- 
icant to lower redshift where they may interfere with attempts to 
measure the power spectrum of ionization fluctuations. The reverse 
is also true. 

6.2 Spectral dependence 

We next imagine using the temperature fluctuations to constrain 
the X-ray source spectra. This should affect the temperature fluctu- 
ations on intermediate scales, where heating fluctuations dominate. 
Increasing the hardness of the spectrum increases the fraction of 
more energetic photons, which have longer mean free paths. This 
should further smooth the temperature fluctuations and suppress 
power on small scales. 
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Figure |9] shows the power spectra at 2: = 15 (chosen to 
maximise the distinctive features of the temperature fluctuations) 
for source spectra as ~ 1.5 (mini-quasars), as = 1.0 (SNR), 
and as ~ 0.5. We see that the spectra alter the most on scales 
A: « 0.1 — 10 Mpc"^. The two main signatures are the change in 
amplitude and shift in the position of the trough. Both of these oc- 
cur because increasing the slope of the spectrum, with fixed total 
luminosity, increases the number of soft X-rays and so increases 
the heating in smaller scales. The trough (or sign change in P^2 ) 
shifts by Ak ~ 2 Mpc~^ for Aq — 0.5, an effect that might be 
observable were it not located on small scales k « 5Mpc~^. The 
amplitude change at the peak is more observable but is also degen- 
erate with modifications in the thermal history, making this a very 
challenging measurement to perform in practice. 

Referring back to our discussion of the time evolution of the 
signal, we see that this sort of variation is similar to the effect 
of changing the thermal history. However, the exact shape of the 
spectrum is determined by the form of gr, and hence Wx- These 
do encode distinct information about the source spectrum. Conse- 
quently, precision measurements of the 21 cm power spectrum at 
high z could constrain the X-ray source spectrum. 

We can also seek to constrain the X-ray spectrum by looking 
at the regime where fluctuations in the Lya flux dominate the 21 
cm signal. The inclusion of Lya photons generated by X-ray ex- 
citation of HI (in addition to those redshifting into the Lyman res- 
onances) modifies the shape of the power spectrum significantly. 
This is easy to see by referring back to Figure [3] There we plotted 
Wa,->c{k), for the case of stellar emission, and Wx{k), which de- 
termines the fluctuations in the X-ray flux. If we allow both stars 
and LyQ photons produced from X-rays to contribute to the Lya 
flux, then the resulting spectrum of fluctuations is determined by a 
weighted c ombination of these Wg (k) and Wx (k). In our model, 
as in that of IChen & Miralda-Escudd ilOO& i. the Lya flux is domi- 
nated on small scales by the X-ray contribution and on large scales 
by the stellar contribution. Thus the resulting weighting function 
most closely resembles Wx{k) with significant power on small 
scales. 

Figure[TO]shows the effect on the power spectrum at 2: = 20, 
when temperature fluctuations are negligible, of including the dif- 
ferent contributions to the Lya flux. On intermediate scales (fc ~ 
1 Mpc~^) there is clearly significantly more power when X-ray ex- 
citation do minates Lya production compared to stellar production. 
As noted in Chuzh ov et al. 112006'), this provides a means for dis- 
tinguishing between the major source of Lya photons during the 
time of the first sources. We note that the shape of the spectrum 
is somewhat sensitive to the spectral index of the X-ray sources - 
with the variation being similar to between the stellar -1- X-ray and 
X-ray only curves. Thus isolating the 21 cm fluctuations from the 
Lya flux variations could also constrain the X-ray spectrum of the 
first sources. 

6.3 Effects of X-ray background 

We now explore the effect of modifying the X-ray luminosity of 
our sources. We have so far taken /x = 1 in our analysis, but 
constraints on the high redshift X-ray background are weak giving 
us significant freedom to vary fx, which parametrizes the source 
luminosity. As an example, for our model A, values of fx ^ 10'' 
are easily possible without X-ray or coUisional ionization of the 
IGM violating WMAP3 constraints on r at the 2-sigma level. In 
Figure [TT] we show the time evolution of the 21 cm fluctuations 
for model A, taking fx = 0.1, 1, and 10. This serves to illustrate 



0.1 
" 102 





'1 ' ' ' 1 ' 


' ■ ■■■ 












r J 
/ 




Star + 
X-ray 


X— ray 

Star ■ 


:(b) 


'1 1 ' 


1 ' ' ' ' 





































10-3 0.01 0.1 1 10 102 

k [Mpc-'] 



Figure 10. {a): |fi,| Aj^j^ (fc). We consider the following sources of Lya 
emssion: stellar only (solid curve). X-ray excitation only (dotted curve), 
stellar+X-ray excitation (dashed curve). All curves are calculated at ^ = 20 
and have been normalised to the stellar only case, to compensate for differ- 
ent mean values of Xa. We assume X-ray emission from starburst galaxies. 
Also plotted is Ift,] A^^ (thin solid curve), (b): \%\^ t^^ ^{k). Same line 
conventions as in fa). 



the effect of late or early X-ray heating and illustrates the range of 
uncertainty in making predictions. 

Earlier heating (dashed curve) causes the temperature fluctua- 
tions to become important at higher redshift, cutting into the region 
of Lya fluctuation. This will make the 21 cm signal more compli- 
cated as temperature and Lya fluctuations contribute over a similar 
range of redshifts. However, early heating also means that temper- 
ature fluctuations become unimportant for the 21 cm signal at late 
times improv ing the prospects for extracting cosmolo gy from the 
21 cm signal jMcOuinn et alj2005lJSantos et al.l2005l) . In contrast, 
late heating (dotted curve) allows a clearer separation between tem- 
perature and Lya fluctuations, but means temperature fluctuations 
are likely to be important during the beginning of reionization. This 
will complicate the extraction of information about HII regions as 
reionization gets underway. 

Clearly there is considerable uncertainty as to the behaviour 
of the 21 cm signal at high redshifts due to our poor understanding 
of the source populations. Viewed another way, measurement of 
the evolution of the 21 cm signal could provide useful constraints 
on the X-ray background at high redshift. This is important as ef- 
forts to observe the diffuse X-ray background are complicated by 
technical issues of calibration. We also note that for weaker X-ray 
heating other sources of heating, especially shock heating, may be- 
come important. 

Finally we remind the reader that our model is applicable in 
the IGM outside of ionized HII regions. If heating occurs late, so 
that temperature fluctuations are important as HII regions become 
large, then it will be important to extend this model if accurate pre- 
dictions of the 21 cm signal during reionization are to be made. It 
will also be important to include these temperature fluctuations into 
simulated predictions of the 21 cm signal. 
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Figure 11. (a): Redshift evolution of ITblArJfc) at k = 0.1 Mpc"! for 
model A, but with fx = 0.1 (dotted curve), 1.0 (solid curve), and 10 
(dashed curve) (b): Redshift evolution of |T(j| A^2 (fc). Same Une conven- 
tions as in (a). 



Figure 12. Predicted l-cr errors on ITjI Ay^^ (k) for an SKA like instrument 
(see text for details). We compare spectra for model A (thick curves) and 
model B (thin curves). Modes with k < ^foreground (shown by a vertical 
dashed line) will probably be lost during foreground cleaning, (a): z = 13. 
(b): z = 15. 



7 OBSERVATIONAL PROSPECTS 



We now turn to the important question of observing the fea- 
tures outlined above. The first generation of 21 cm experiments 
(PAST, LOFAR, MWA) will be optimised to look for the sig- 
nature of HII regions at redshifts z < 12. Their sensitivity 
decreases rapidly at redshifts z > 10 jMcOuinn et al. I I2OO5I : 
iBowman. Morales & Hewit3l2006h and so they are unlikely to be 
able to detect the effects of inhomogeneous heating. The proposed 
successor to these instruments, the SKA, is still under design, but 
its fiducial specifications should allow the 2: > 12 regime to be 
probed. In this section, we will consider using an SKA type exper- 
iment to observe 21 cm fluctuations at z = 13 and z — 15 and 
calculate the achievable precision. 

Before this, we must make the necessary caveats concern- 
ing foregrounds. Foregrounds for 21 cm observations include 
terrestrial radio interference (RFI), galactic synch rotron emis- 
sion, radio recombination lines, and ma ny others dOh & Mack] 
20031; iDi Matteo, Ciardi & MiniatH |2004 §9 of iFurlanetto etd 



20061) . Typical foregrounds produce system temperatures Tsj 



1000 K, compared to a signal measured in mK. These fore- 
grounds increase rapidly as we move to lower frequency, mak- 
ing their removal an even greater concern for high redshift ob- 
servations than low ones. Although techniques for foreground re- 
moval are well grounded, their effectiveness has yet to be tested. 
In the analysis that follows, we assume that foreground removal 
can be effected by exploiting the smoothness of foregrounds 
in frequency space Iz aldarriaga, Furlanetto & Hemguisll 20041: 
Santos, Coorav & Knox 2005; Morales & Hewitll |2004 iMoralej 



2005; McOuinn et al. 2005; Wang & Hu 2006). 

Figure [12] shows predicted 1 — a error bars on Atj, (fc) at 
z = 13 and z — 15 for model A and B. We assume an SKA like in- 
strument with a total effective area Atot = 1 km^ distributed over 
5000 antennae in a 5 km core, bandwidth 5=12 MHz, minimum 
baseline Dmin ~ 10 m, and an integration time tint = 1000 hr. 



bins of width Ak = k/2. We assume that foregrounds can be re- 
moved exactly, but that this also removes cosmological informa- 
tion on scales exceeding the bandwidth of the observations, so that 
modes with k < feforcgro und ~ 0.025 Mpc~^ (indicated by verti- 
cal dashed lines) are lost jMcOuinn et al 

With these caveats, observations could measure At^ (k) ac- 
curately over the range k ~ 0.025 — 2 Mpc~ . The precision is 
more than adequate to distinguish between model A and B. De- 
tecting the characteristic peak-trough signature of St is difficult, as 
the trough typically occurs on small scales where the uncertainty 
is large. However, it should be possible to detect the peak and the 
beginning of the decline. We note that detection of the trough is 
necessary to unambiguously determine /3t < and so show that 
Ts < T-y. Without this it is not simple to distinguish between the 
two cases exemplified by the z — 14 curve, which has no trough, 
and the z = 15 curve, which does, in Figure [8] No similar con- 
fusion occurs when the reduction in power caused by /3t < is 
obvious, as in the z — 15 model B case. 

From the point of view of constraining the spectra of X-ray 
sources, the precision is adequate for distinguishing between the 
different curves of Figure [lO] Whether the effect of the spectrum 
can be separated out from different thermal histories is an open 
question, which deserves future study. 

Throughout this work, we have ignored the effect of the HII 
regions on the 21 cm power spectrum. While this is reasonable at 
high redshifts, this approximation will begin to break down as the 
filling fraction of ioniz ed regions increases. The bubble model of 
IFurlanetto et al. I ( I2OO4I) predicts that these bubbles remain at sub- 
Mpc sizes while Xi < 0.1. Consequently, we naively expect con- 
tamination of the signal by these bubbles to be confined to small 
scale modes with k > 1 Mpc~^ that will be very difficult to detect. 
Exploring the detailed interaction between temperature and neutral 
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fraction fluctuations is beyond the scope of this paper, but may be 
important for detailed predictions of the 21 cm signal at the begin- 
ning ofreionization. 

ISantos &Cooravl ( l200dh have considered the extraction of as- 
trophysical and cosmological parameters from 21 cm observations 
in the period of the first sources. They assumed that gas tempera- 
ture fluctuations showed no scale dependence gT{k,z) = griz) 
and argued that extracting astrophysical information using an SKA 
like instrument is difficult but feasible. We expect the scale depen- 
dent temperature fluctuations that we have investigated to both help 
and hinder parameter estimation. Figure[T2lshows that it should be 
possible to resolve individual features imprinted in the power spec- 
trum by temperature fluctuations. These features provide additional 
leverage in extracting astrophysical parameters. However, the shape 
of the power spectrum evolves rapidly in our model, making bin- 
ning of different redshift data more difficult. 



8 CONCLUSIONS 

X-ray production by an early generation of stellar remnants is 
widely regarded as the most likely candidate for heating the IGM 
above the CMB temperature from its cool adiabatic level. This heat- 
ing has often been treated as uniform, as the mean free path of hard 
X-rays in the early Universe is comparable to the Hubble scale. 
We ha ve relaxed this assumptio n and, by expanding on the formal- 
ism of lBarkana & Loebl j2005bh , calculated the temperature fluctu- 
ations that arise from the inhomogeneous heating. The spectrum of 
fluctuations in Tk is significantly larger than that predicted from 
uniform heating, peaking on scales k ~ 0.1 Mpc~^. This allowed 
us to examine the redshift range about Zh, where Tk = T7, and 
show that there is a window of width Az ~ 1 in which the IGM 
will contain pockets of gas both hotter and colder than the CMB. 
This has implications for the 21 cm signal, which will be seen in a 
mixture of absorption and emission within this window. 

The best hope for observing the temperature evolution be- 
fore reionization is through 21 cm observations of neutral hydro- 
gen. Systematic effects arising from foregrounds are like ly to pre- 
vent interferometers from measuring TJ, directly (Furlanetto et al.l 
althoug h several alternative methods for obtain ing Tt have 
been proposed tearkana & Loebll2005al :l Cooravir2006[) . Thus care- 
ful analysis of brightness fluctuations will be required to extract 
astrophysical information. Fluctuations in Tk lead to fluctuations 
in Tb, which contain information about the thermal history and the 
nature of the heating sources. We have calculated the 21 cm power 
spectrum arising from inhomogeneous X-ray heating and shown 
that it has considerable structure. In the regime where gas tem- 
perature and Lya flux fluctuations compete, we expect a trough- 
peak-trough structure in AT^{k). Once Tk fluctuations dominate, 
but while Tk < T^,, we see a peak-trough structure. As the gas 
heats, this structure is lost as the trough moves to unobservable 
small scales while the peak decreases and finally vanishes once 
Tk ^ T^- Extracting astrophysical information cleanly will be 
challenging, but the information is there. 

It is important to notice that the difference between uniform 
and inhomogeneous heating is large. Observations with the SKA 
should be able to distinguish these two cases and indicate whether 
X-ray heating is important. If it is possible to perform an angular 
separation of Ptj, then observing P^2 < is a clear indicator 
that Tk < T-,. Ideally, one would extract the quantity /3t, but this 
requires fitting of other parameters and so is a less direct (but more 
conclusive) observational feature. 



Additionally, the spectra of the X-ray sources imprint infor- 
mation on the Tk fluctuations. This may be observed in the 21 
cm power spectra, where it shifts the critical scale at which P^2 
changes sign, or during the regime in which Lya fluctuations dom- 
inate, where it modifies the shape of the power spectrum. The tem- 
perature fluctuations that we have calculated lead to a 21 cm signal 
that extends down to relatively low redshifts. This opens an oppor- 
tunity for future 21 cm radio arrays to probe the thermal history 
prior to reionization. Including temperature fluctuations makes the 
21 cm signal significantly more complex, adding information, but 
further raises the question of how best to separate out that informa- 
tion. 

In this paper, we have ignored the contribution from Poisso n 
fluctuations in the source distribution tearkana & Loebl l2005bl) . 
While calculating it requires only a straightforward extension of 
the Barkana & Loeb (2005b) formalism, performing the time in- 
tegrals necessary to convert heating fluctuations into temperature 
fluctuations is non-trivial. We have estimated the amplitude of these 
Poisson temperature fluctuations and find them (in our models) to 
be subdominant at all redshifts. This is largely because there are 
many more sources at the lower redshifts where temperature fluc- 
tuations are important. In theory, high-precision 21 cm observations 
can separate these Poisson fluctuations from fluctuations correlated 
with the density field. The Poisson contribution could then be used 
to probe the distribution of sources, for example, by distinguish- 
ing between highly biased mini-quasars and less biased star-burst 
galaxies, producing the same global X-ray luminosity. 

In our analysis we have taken /x = 1, corresponding to nor- 
malising the X-ray luminosity per unit star formation to that ob- 
served in the local universe. In truth, this assumption is highly spec- 
ulative and the value for fx is extremely uncertain. We have inves- 
tigated the effects of changing fx and find that it alters the details 
of the thermal evolution significantly. Taking fx ~ 0.1, for ex- 
ample, shifts the point where 21 cm brightness fluctuations change 
from being dominated by Lya fluctuations to gas temperature fluc- 
tuations from z ~ 17 to z « 15. Setting fx = 10 increases the 
redshift of this transition to 2; ~ 19. For values of fx ^ 0.1, 
we find a clear separation between 21 cm brightness fluctuations 
sourced by gas temperature and Lya fluctuations. Increasing fx 
also increases the redshift at which Tk ^ T^,, so that gas tem- 
perature fluctuations become irrelevant for the 21 cm signal. Addi- 
tionally, small values for fx will increase the contribution of other 
heating mechanisms such as shock heating. All of this suggests that 
measuring the time evolution of the 21 cm signal (as in Figure [TT] 
for example) would enable fx to be constrained. Unfortunately, 
until these observations are made it is difficult to predict the ther- 
mal history before reionization with any certainty. 

We have shown that the 21 cm signal at high z will con- 
tain significantly more structure than has previously been con- 
sidered. Temperature fluctuations produce an interesting interplay 
with other sources of 21 cm anisotropy as /3t < when Tk < T-,. 
Furthermore, for reasonable heating scenarios, the effect of temper- 
ature fluctuations persist well into the regime that will be probed by 
second generation low-frequency arrays, such as the SKA. Thus, 
prospects for probing the thermal history before reionization via 
observations of the redshifted 21 cm line seem promising. 
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